New mass limit for white dwarfs: super-Chandrasekhar type la supernova as a new 

standard candle 



43 



CO 
C3 



> 
in 



is 



iMCS: 97.60.Bw, 97.20.Rp, 98.80.Es, 97.10.Ld, 71.70.Di 



Upasana Das, Banibrata Mukhopadhyay 
Department of Physics, Indian Institute of Science, Bangalore 560012, India 
upasana@physics.iisc.ernet.in , bm@physics.iisc.ernet.in 

Type la supernovae, sparked off by exploding white dwarfs of mass close to Chandrasekhar 
limit, play the key role to understand the expansion rate of universe. However, recent obser- 
vations of several peculiar type la supernovae argue for its progenitor mass to be significantly 
super-Chandrasekhar. We show that strongly magnetized white dwarfs not only can violate the 
£f~^ • Chandrasekhar mass limit significantly, but exhibit a different mass limit. We establish from foun- 

dational level that the generic mass limit of white dwarfs is 2.58 solar mass. This explains the origin 
f"^ , of over-luminous peculiar type la supernovae. Our finding further argues for a possible second stan- 

£SJ ' dard candle, which has many far reaching implications, including a possible reconsideration of the 

| ' expansion history of the universe. 

I—* 

in 

Introduction.— Recently, some peculiar type la supernovae: e.g. SN 2006gz, SN 2007if, SN 2009dc, SN 2003fg, have 
fV^ been observed [l| with exceptionally higher luminosities but lower kinetic energies. The kinetic energy is proportional 
' to the difference between the obtained nuclear energy which arises from the synthesis of elements in the explosion 
through fusion and the binding energy of the white dwarf. Most of the light-curves of above mentioned peculiar 
supernovae appear over-luminous, slow-rising, which do not allow them to be calibrated as standard candles. This 
questions the use of all type la supernovae in measuring distances of far away regions and hence unraveling the 
expansion history of the universe. However, assuming the progenitor to be a highly super-Chandrasekhar mass white 
dwarf reproduces the low kinetic energies and thus velocities seen in the above supernovae Q . This is because a larger 
mass implies a larger binding energy of the star and hence a smaller velocity for the same and/or higher luminosity 
(due to nuclear fusions) than that observed in a standard type la supernova. The progenitor masses required to 
explain the above supernovae lie in the range 2.1 — 2.8M Q , when M Q being the mass of sun, subject to the model 
chosen to estimate the nickel mass [3-[B] ■ Now naturally the following vital questions arise. Is there any fundamental 
basis behind the formation of a highly super-Chandrasekhar white dwarf? How to address the significant violation 
of Chandrasekhar mass limit? What is the ultimate mass limit of a white dwarf? Here we plan to address all the 
Q\ ' above issues by exploiting the effects of magnetic field in compact objects. This will lead, as we will show, to a 
I/"") , natural explanation of the so called "peculiar" type la supernovae. This might eventually lead these supernovae to 
be considered as altogether new standard candles. This has many far reaching implications, including a possible 
, reconsideration of the expansion history of the universe, 
(f) ■ Before proceeding further, let us recall the physics of a white dwarf and its link to the type la supernova. When a 
t-H | star exhausts its nuclear fuel, it converts to either of the three compact objects: white dwarf, neutron star or black 
hole, depending on the initial mass of the evolving star. It is generally believed that the fate is a white dwarf when the 
mass of the initial star in its main sequence is Mms~5Mq. Such a main sequence star undergoing collapse leading to 
a small volume consists of a lot of electrons. Being in a small volume many such electrons tend to occupy same energy 
states, hence making them to degenerate electrons, as the energy of a particle depends on its momentum which is 
determined by the total volume of the system. However, being a fermion an electron obeys Pauli's exclusion principle 
which says that no two fermions can occupy the same quantum state. Hence, once up to the Fermi level, which is the 
maximum allowed energy of a fermion, is filled by the electrons, there is no available space for the remaining electrons 
in a small volume of a collapsing star, which expels the electrons to move out leading to an outward pressure. In a 
white dwarf, the inward gravitational force is balanced by the force due to outward pressure created by to degenerate 
electrons. Chandrasekhar in one of his celebrated papers Q showed that the mass of a white dwarf cannot be more 
than 1.44Mq which sets the famous Chandrasekhar mass limit of a white dwarf. 

If a white dwarf having mass close to the Chandrasekhar limit gains more mass (e.g. by accretion, when the mass 
is supplied by a companion star of the white dwarf), then its mass exceeds the Chandrasekhar limit, which leads to 
a gravitational force stronger than the outward force that arises due to the degenerate electrons. Hence, this leads to 
the contraction of the white dwarf and a subsequent increase of its temperature, which is favorable for the initiation 
of fusion reactions again. If the white dwarf mostly consists of carbon and oxygen, namely the carbon-oxygen white 
dwarf (which is commonly the case), then nuclear fusion of carbon (and oxygen) takes place. Subsequently, within 
a few seconds, a substantial fraction of the white dwarf matter undergoes a runaway reaction which releases huge 
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energy ~ 10 51 erg to unbind it in an explosion, namely type la supernova explosion. This eventually leads to a complete 
gravitational collapse of the star without leaving any remnant. 

As all the commonly observed type la supcrnovae are produced by (almost) the same mechanism, namely the mass 
of the progenitor white dwarf exceeding Chandrasekhar limit and subsequent processes, the underlying variations 
of luminosity as functions of time, namely light-curves, appear alike for all the explosions. All of these supernovae 
exhibit consistent peak luminosity, the relation between the peak luminosity and width of the light-curve [1, Hj] , due to 
the uniform mass of white dwarfs (Chandrasekhar limit) which finally explode, e.g., because of the accretion process. 
Very importantly, the apparent stability of this value helps the underlying supernovae to be used as standard candles 
in order to measure the distances to their host galaxies. Since these supernovae are exceptionally bright, they can be 
observed across huge cosmic distances. The variation of their brightness with distance (or redshift) is an extremely 
important tool for measuring various cosmological parameters, which in turn shed light on the expansion history of 
the Universe. Their enormous importance is self evident and was brought into the prime focus by the awarding of the 
Nobel Prize in Physics in 2011, for the discovery (made possible by the observation of distant type la supernovae) 
that the universe is undergoing an accelerated expansion [1(1 HH • 

Now we move on to our goal of establishing a new mass limit for super-Chandrasekhar white dwarfs, of whose 
formation there is no understanding from foundational level — a caveat behind the hypothesis of super-Chandrasekhar 
progenitor for the peculiar type la supernovae, raised by the earlier authors 0. In fact they emphasized on the 
pursuit of theoretical studies in order to assess the hypothesis. Although, based on a numerical code for stellar binary 
evolution, the rotating white dwarfs are suggested to hold mass up to 2.7M Q [12], a foundational level calculation is 
missing and there is no estimate of mass limit of such a star either. 

In this letter, we show that (highly) magnetized white dwarfs not only can have mass ~ 2.6M Q , but also exhibit its 
ultimate limit of mass. Hence, we propose a fundamentally new mass limit for white dwarfs, which eventually helps 
in explaining light-curves of peculiar type la supernovae. This may further lead to establishing these supcrnovae as 
modified standard candles for distance measurement. 

The motivation behind our approach lies in the discovery of several isolated magnetized white dwarfs through the 
Sloan Digital Sky Survey (SDSS) with surface fields 10 5 — 10 9 G [HI, [l4|. Hence their expected central fields could 
be 2 — 3 orders of magnitude higher. Moreover, about 25% of accreting white dwarfs, namely cataclysmic variables 
(CVs), are found to have high magnetic fields 10 7 — 10 8 G [HI]. 

Equation of state.— As the starting choice is the magnetized white dwarf, we first recall degenerate electrons 
under the influence of magnetic field which are known to be Landau quantized Larger the magnetic field, smaller 
is the number of Landau levels occupied [ItJ (see supplement [U|)- Recent works H(| establish that Landau 

quantization due to strong magnetic field modifies the equation of state (EoS) of the electron degenerate gas, which 
results in a significant modification of the the mass-radius relation of the underlying white dwarf. Interestingly, these 
white dwarfs are found to have super-Chandrasekhar masses. The main aim of this letter is to obtain the maximum 
possible mass of such a white dwarf (which is magnetized), and therefore a (new) mass limit. Hence we look for the 
regime of high density of electron degenerate gas and the corresponding EoS, which further corresponds to the high 
Fermi energy (Ep) of the system. This is because high density corresponds to high momentum, which implies high 
energy (see supplement [HI]). Note that the maximum Fermi energy (Ep max ) corresponds to the maximum central 
density of the star. Consequently, conservation of magnetic flux (technically speaking flux freezing theorem which is 
generally applicable for a compact star) argues for the maximum possible field of the system, which implies that only 
the ground Landau level will be occupied by the electrons. For the expressions of density, pressure and the EoS for 
such a highly magnetized system, see supplement [l8j . Hence, in the limit of Ep >> m e c 2 , when m e is the mass of 
the electrons and c the speed of light, for a given magnetic field exhibiting the system of one Landau level, the EoS is 

P = K rnP 2 , (1) 

when P and p are respectively the pressure and density of the gas and the constant K m is given by 

_ m e c 2 ir 2 \ 3 e 
{li e m H yB D 

where A e = h/m e c, the Compton wavelength of electron, h the Planck's constant h divided by 2-7T, p e the mean 
molecular weight per electron, m# the mass of hydrogen atom and Bp, the magnetic field in the units of 4.414 x 10 13 G. 

Mass limit of white dwarfs.— Now following the Lane-Emdcn formalism [4J, we obtain the mass of the magnetized 
white dwarf 
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and the corresponding radius 



when p c is the central density of the white dwarf supplied as a boundary condition in addition to the condition that 
dp/dr = at r = 0, when r is the radial distance from the center of the star such that at the surface r = R, and G is 
the Newton's gravitation constant. See supplementary information [l8| for detailed calculations. 

Now the expression of p c for a one Landau level system in the limit Ep = Ep max » m e c 2 is given by (see 
supplement [18j |) 

^emH „3/2 ( r-s 

Substituting p c from equation ([5]) in equation (j3|), we obtain the mass of the white dwarf, independent of p c and Bd, 
given by 

/ hc\ 3/2 1 10.312 w 

\2GJ {p e m H y n 2 e 

and from equations ([2j), (jlj and (0 we obtain the radius 

( * V3 he \ 1/2 

R= { ^ (p emH )^ G ) AT 7 ->0bs Aj -kx>, (?) 

which set the new limits for mass and radius. Note that the Chandrasckhar limit also corresponds to p — > oo and 
R — > 0. For p e = 2 which is the case of a carbon-oxygen white dwarf 

M w 2.58M Q . (8) 

Interestingly, while high magnetic field introduces anisotropic effects into the star tending it to be an oblate spheroid, 
this does not affect the limiting mass as the corresponding radius tends to zero. However, for lighter white dwarfs 
with finite radii, super- Chandrasckhar mass would have been achieved at a lower field, if the star is appropriately set 
to be a spheroid rather than a sphere, as justified earlier pjl f20| . 

Scaling behaviors of mass and radius of the white dwarf with its central density.— Now we provide general scaling 
laws for the mass and radius of the white dwarf describing their variations with its central density (as is known for 
non-magnetized white dwarfs proposed by Chandrasckhar) and magnetic field strength. By Lane-Emden formalism 

MocK 3/2 p7^, (9) 

where K depends on Bo and p c for a magnetized white dwarf (see supplement [Hj]). For the extremely high density 
regime of the white dwarf, the EoS reduces to a polytropic form P = Kp T (can be verified from supplement [IH), 
when r = 1 + 1/n, is the polytropic index. In this regime, r = 2, consequently n — 1 and K — K m oc B^ 1 cx p c 2 ^ 3 , as 
shown by equations (|I.5|) . ([2]), <j5j) . This finally renders M in equation (|9]) to be independent of p c , as already shown 
by equation (O. Similarly, the scaling of radius is obtained as 

RocK 1 /^^, (10) 

which reveals R cx p c for n = 1, as already shown by equation ([?])• 

Justification of high magnetic field in white dwarfs.— So far we have employed Landau quantization in search of a 
new mass limit giving rise to super-Chandrasekhar white dwarfs. However, the effect of Landau quantization becomes 
significant only at a high field ~ Bp> x 10 13 G = B cr . How can we justify such a high field in a white dwarf? 

Let us consider the commonly observed phenomenon of a magnetized white dwarf accreting mass from its companion. 
Now the surface field of an accreting white dwarf, as observed, could be ~ 10 9 G << B cr [13l |. Its central field however 
can be several orders of magnitude higher ~ 10 12 G, which is also less than B cr . Naturally, such a magnetized CV, 
commonly known as a polar, still lies on the mass radius relation obtained by Chandrasekhar. However, in contrast 
with Chandrasckhar's work (which did not include magnetic field in the calculations), we will see that a non-zero 
initial field in the white dwarf, however ineffective for rendering Landau quantization effects, will prove to be crucial in 
supporting the additional mass accumulated due to accretion. As the magnetized white dwarf accretes mass, its total 
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FIG. 1: Equations of state — the solid line represents Chandrasekhar's equation of state (corresponding to zero magnetic 
field and hence infinitely many Landau levels), the dotted and dashed lines represent the 5- level (corresponding to very 
strong magnetic field) and f-level (corresponding to ultimate equation of state for extreme magnetic field) systems of Landau 
quantization respectively. EFmax = 200m e c 2 . 



mass increases which in turn increases the gravitational power and hence the white dwarf contracts in size due to the 
increased gravitational pull. However, the total magnetic flux is conserved in such a process and, hence, as a result 
of the above decrease in size of the star, the central (as well as surface) magnetic field also increases. Here we are 
interested in the evolution of the central field, since it is this field which is primarily responsible for rendering super- 
Chandrasekhar mass to the white dwarf, as justified in Q. Since accretion is a continuous process, the deposition 
of matter on the surface of the white dwarf, followed by its contraction and subsequent increase of magnetic field, 
continues in a cycle. In such a process, eventually the central magnetic field could exceed B cr . As a result, the 
EoS of the electron degenerate matter gets modified as shown in Figure 1. Hence, the inward gravitational force is 
balanced by the outward force due to this modified pressure and a quasi-equilibrium state is attained. In this way, a 
very high magnetic field is generated, which in turn prevents the white dwarf from collapsing, thus making it more 
massive. Subsequently, with the continuation of accretion the white dwarf approaches the new mass limit ~ 2.58Mq, 
as obtained above, which sparks off a violent thermonuclear reaction with further accretion, thus exploding it and 
giving rise to a super- Chandrasekhar type la supernova. The evolution of the mass-radius relation of a polar into 
that of a super-Chandrasekhar white dwarf of the maximum possible mass is shown in Figure 2, along with a few 
typical mass-radius relations for different magnetic field strengths describing possible stars in intermediate equilibrium 
states. The ultimate white dwarf, corresponding to the maximum mass ~ 2.58Mq, lies on the mass-radius relation for 
a one Landau level system, but the intermediate white dwarfs having weaker magnetic fields correspond to multilevel 
systems. The one Landau level system corresponds to the central magnetic field 8.8 x 10 17 G. The intermediate 
systems of 200-level and 50124-level correspond to central magnetic fields 4.4 x 10 15 G and 1.7 x 10 13 G respectively 
(see supplement (l8T | for relevant formula), when Ep max = 200ra e c 2 . This value of Ep max is found to produce the 
theoretical mass limit with good numerical accuracy. Note that one can in principle go upto higher Ep max , which will 
lead to a further decrease in the radius of the white dwarf keeping the mass practically same. In order to construct 
this evolutionary track, we consider the values of p c corresponding to the density at the ground-to-first Landau level 
transition of the respective EoSs, since only this choice leads to the maximum possible mass. We emphasize here that 
the range of masses for the super-Chandrasekhar progenitors obtained from observations is not very strict. Hence, 
if the new mass limit obtained by us is taken into account, one could possibly do away with the mass distribution 
altogether. However, if the distribution is indeed real, then it could be attributed to the difference in accretion rates 
found in different CVs. 

We thank A.R. Rao of TIFR for insightful suggestions. B.M. acknowledges partial support through research Grant 
No. ISRO/RES/2/367/10-11. U.D. thanks CSIR, India for financial support. 
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FIG. 2: Mass-Radius relations — the pure solid line represents Chandrasekhar's result and the one marked with filled circles 
represents the evolutionary track of the white dwarf with the increase of magnetic field. The dot-dashed, dotted and dashed lines 
represent the white dwarfs with 50124-level, 200-level and 1-level systems of Landau quantization respectively (corresponding 
to increasing magnetic fields). Eir max = 200m e c 2 . 
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Supplementary information for the letter 
I. BASIC EQUATIONS FOR ONE LANDAU LEVEL SYSTEM OF DEGENERATE ELECTRONS 

The Fermi energy level (Ep) of a Landau quantized system is given by 

E% = PF {vfc 2 +mlc i {l + 2 V B D ), (1.1) 

when pf{v) is the Fermi momentum of vth Landau level (y — 0, 1,2 ), m e the mass of electrons, c the speed of 

light, Bp, the magnetic field in units of 4.414 x 10 13 G. From the condition that pf(v) 2 > 0, the maximum number of 
occupied Landau levels is given by [H 



/ -Efmax V 

\ m e c 2 j 



2B D 



(1.2) 



For one Landau level system, when only ground Landau level (y = 0) is occupied, v m = 1. Similarly, for two level 
system, when ground (v = 0) and first (y = 1) levels are occupied, v m = 2, and so on. Hence, for one Landau level 
system, density (p) and pressure (P) of the electron degenerate gas 0] are given by 



p e m H B D If Ep \ 2 p e m H 



p = ^aTV J 1 2M! Bc£f for >:> meC ' (L3) 



B D m e c 2 E F ^/E F -mlc* \ E F + y/E*. - m'c* \\ 
4^ 2 A 3 ^ (m e c 2 ) 2 \ m e c 2 J J ' 1 ' 

when fi e is the mean molecular weight, m# the mass of proton, A e the Compton wavelength of electron. Now 
eliminating Ep from equations (|I.3[) and p.4|) for an arbitrary Ep, we arrive at the equation of state (EoS) for a one 
level system given by 



when Q = p e mHB D /2ir 2 \^,. 

Following previous work 0] , we now approximate the above EoS and EoSs for any other v m by a polytropic relation 
P = Kp T such that the polytropic index T = 1 + 1/nis piecewise constant in different density ranges and K being a 
dimensional constant. This will prove to be useful in order to understand scaling behaviors of mass and radius of the 
white dwarf with its central density. 

Now in terms of Q, equation (|I.3[) can be rewritten as 

P = Q^- (L6) 
In the limit Ep >> m e c 2 , we note that, p » Q. Thus in this limit, equation (|I.5[) reduces to 

/J 2_ g 2 ln r^iy (L7) 



2Qfi e m H 

The logarithmic term is much smaller than the first term in the above equation and hence by neglecting it we obtain 

P=^^P\ (L8) 
which corresponds to the polytropic EoS with V = 2. 
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II. LANE-EMDEN EQUATION AND EXPRESSIONS FOR MASS AND RADIUS 



The underlying white dwarf following above EoS obeys the magnetostatic equilibrium condition 

1 d ( B 2 \ B-VB , s 

- p d-r{ P+ 8;)= F ° + —' ^ 

when r is the radial distance from the center of white dwarf, B the magnetic field in G, B 2 = B-B, F g the gravitational 
force. This equation is supplemented by the estimate of mass (M) within any r given by 

= 4nr 2 p II.2) 

dr 

approximating the star to be spherical. The white dwarf in the present context can be considered in the Newtonian 
framework and the magnetic field therein is nearly constant in the regime of smaller radii (as justified previously 0] ) ■ 
Moreover, at a very large density, as in the limiting case to be considered here, the star becomes so compact as if the 
magnetic field remains constant throughout. Hence, taking above facts into consideration and combining equations 
(III. II) and pi.2[) . the equilibrium condition may be read at any r as 

LA. ( r l^L\ = -^Gp (II 3) 

r 2 dr \ p dr J 

where G is the Newton's gravitation constant. We now briefly recall the Lane-Emden formalism (see, e.g., Q). Let 
us define 

P = PcO n , (11 A) 

where p c is the central density of the white dwarf and 9 is a dimensionless variable and 

r = at, (11.5) 
where £ is another dimensionless variable and constant a carries the dimension of length defined as 

1/2 

(n + l)Kp c n 

a 



(11.6) 



4ttG 

Thus using equations pi.4[) , (|II.5|) and (|II.6|) along with the polytropic form of EoS, equation (III. 31) reduces to 

H ( eJ i) - (IL7) 

which is the famous Lane-Emden equation. Equation pi.7[) can be solved for a given n, along with the boundary 
conditions 

0(£ = 0) = 1 (II.8) 

and 

'd6 s 



di 



= 0. (II.9) 

Note that for n < 5, 8 becomes zero for a finite value of £, say £i, which basically corresponds to the surface of the 
white dwarf such that its radius 



i? = a£i. (11.10) 
Also by combining equations pl.2j) . pl.4|) . (|II.5|) and (|II.7|) we obtain the mass of the white dwarf 
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M = 4ira 3 Pc I £ 2 <T df. (11.11) 



n 
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